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Low dimensionality, broken symmetry and easily-modulated carrier concentrations provoke novel 
electronic phase emergence at oxide interfaces. However, the spatial extent of such reconstructions - i.e. the 
interfacial "depth" - remains unclear. Examining LaAlOs/SrTiOs heterostructures at previously unexplored 
carrier densities — 6.9 X 10^^ cm"^, we observe a Shubnikov-de Haas effect for small in-plane fields, 
characteristic of an anisotropic 3D Fermi surface with preferential d^^^yz orbital occupancy extending over at 
least 100 nm perpendicular to the interface. Quantum oscillations from the 3D Fermi surface of bulk doped 
SrTiOs emerge simultaneously at higher ii2D- We distinguish three areas in doped perovskite 
heterostructures: narrow (<20 nm) 2D interfaces housing superconductivity and/or other emergent 
phases, electronically isotropic regions far (> 120 nm) from the interface and new intermediate zones where 
interfacial proximity renormalises the electronic structure relative to the bulk. 



Ever since the discovery of a conducting channel in LaAlOs/SrTiOs^ and the subsequent observations of 
magnetism^ and superconductivity^ the vast majority of oxide interface research has focussed on synthe- 
sising intrinsically- doped heterostructures featuring narrow conducting channels (<20 nm) with two- 
dimensional carrier densities in the 10^^-10^^ cm"^ range^"^. At such interfaces, it has been shown^'^^"^^ that 
symmetry-lowering and quantum confinement lift the Ti degeneracy, so that the d^y orbital lies at lower energy 
than the dxz,yz orbitals. X-ray absorption spectroscopy^ reveals a band splitting of —50 meV for ~ 10^^ cm"^ 
and theoretical approaches indicate that this increases with n2D> reaching —0.25 eV at 3 X 10^^ cm"^ 
Regardless of the total /12Dj the splitting should gradually vanish below the interface, until the electronic structure 
resembles that of bulk SrTi03 with degenerate dxy,xz,yz orbitals creating a Fermi surface at the centre of the Brillouin 
zone^^ The lengthscale over which this degeneracy is regained - i.e. the total distance over which the interface induces 
electronic reconstruction - remains unknown, despite being a vital prerequisite for building layered 3D oxide devices. 

Probing this lengthscale requires the synthesis of LaAlOs/SrTiOs heterostructures with significantly more 
carriers (and correspondingly deeper conducting channels) than the norm. Previously, high n2D heterostructures 
have only been grown in reducing environments^^'^^, creating bulk-like conducting layers hundreds of microns 
thick {n2D — 10^^ cm"^) in which the broken symmetry of the interface plays no role. However, interfaces with 
5x10^^ < n2D < 5x10^^ cm ~^ have until now remained unexplored: at these intermediate n2D^ electrons "spill 
over" from the interface and begin to occupy states lying deeper within the SrTiOs. The principal focus of our 
work is therefore to track the evolution of the electronic structure and its crossover from 2D interfacial to 3D bulk- 
like behaviour within this range of carrier densities. For n2D — 6.9 X 10^^ cm"^, we report the first instance of 
Shubnikov-de Haas (SdH) oscillations from an ultra-high mobility electron gas (^h ~ 10^ cm^V^s"^) for small 
magnetic fields parallel to the interface. The absence of such oscillations from the low-field perpendicular 
magneto resistance indicates that these carriers originate from an anisotropic 3D Fermi surface (FS); our first- 
principles calculations of the sub -interfacial electronic structure reveal dominant dxz,yz orbital occupancy, which 
is consistent with our experimental data. Superconductivity remains confined within 20 nm of the interface, 
while the 3D FS characteristic of bulk doped SrTiOs gradually emerges with increasing n2D- Together, our results 
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imply the existence of a region below the interface whose electronic 
structure differs from that of the bulk, with a minimum thickness of 
100 nm imposed by the cyclotron radius. This discovery has import- 
ant implications for oxide devices seeking to functionalise interfacial 
electronic reconstructions. 

Results 

During sample growth, three mechanisms exist for carrier- doping 
the LaAlOs/SrTiOs interface: (a) intrinsic self-doping via the polar 
catastrophe^^ (leading to a maximum — 3.3 X lO^'* cm"^), (b) 
oxygen vacancy doping^^ (contributing 2e~ per O^" vacancy) and (c) 
cation intermixing^ (an unbalanced switching of La^^ for Sr^^ and 
AP^ for Ti^^). Since our principal aim is to explore the evolution of 
the electronic structure for > 5 X 10^^ cm"^ (far beyond the 
upper limit imposed by the polar catastrophe) and cation intermix- 
ing is difficult to control in a pulsed laser deposition (PLD) chamber, 
we use O^" vacancy doping to achieve the high values necessary 
for this project. To this end, we synthesise LaAlOs/SrTiOs hetero- 
structures at an intermediate O2 pressure (10"^ mbar), without any 
post- annealing procedure (further growth and characterisation 
details may be found in the Methods and Supplementary 
Material). The lack of annealing guarantees a high vacancy 
concentration and hence a large n2Dy while the intermediate growth 
pressure ensures that these vacancies do not penetrate far into the 
SrTiOa substrate. Low pressure growth (10"^ mbar) without anneal- 
Ij^gi6,i7 j^^g previously been shown to result in macroscopic substrate 
conduction, with /I2D ^ 5 X 10^^ cm"^; in contrast, our method of 
synthesis consistently yields heterostructures with as -grown Hall 
carrier densities in the 10^^-10^^ cm"^ range, which we will refer to 
as "series B". For comparative purposes, we have also annealed cer- 
tain heterostructures ("series A") at high O2 pressures, yielding 
~ 10^^ cm"^. A-type interfaces are comparable to the majority of 
those previously studied in the literature^'^'^, in which carrier injec- 
tion is dominated by the polar catastrophe. Both series exhibit coex- 
istent superconductivity (SC) and ferromagnetism (FM), a 
comprehensive analysis of which may be found in ref. 20. For quant- 
itative continuity in the present work, we focus on two specific sam- 
ples A and B, with as-grown ^20 = 2.3 X 10'' cm"^ 6.9 X 10'^ cm"^ 
at T = 0.1 K and SC channel thicknesses d = 18 ± 1 nm, 9 ± 1 nm 
respectively. Sample B has a back gate beneath the SrTiOs: /I2D 
increases to 2.4 X 10'^ cm"^ (d = 19 ± 2 nm) at gate voltage Vg = 
350 V. The heterostructure withstands Vg = 500 V with no discern- 
ible leakage current and the substrate capacitance ~ 1 nF is compar- 
able to values measured in annealed LaAlOs/SrTiOs heterostructures 
with lower /i2D^'^ (see Supplementary section 1). Such conditions can 
only be achieved if the bulk of the SrTiOs substrate is insulating: this 
confirms that vacancies have not penetrated deep into the SrTiOa 
and are restricted to the neighbourhood of the interface. 

We probe the electronic structure and FS geometry using SdH 
oscillations in the magnetoresistance (MR) RxxiH) (Fig. la). Two 
magnetic field orientations are principally considered: H ± (001) 
(H^) and H//[110] (H//), where the [100] directions correspond to 
the crystallographic axes of the SrTiOs substrate and [001] points 
out-of-plane. Sample A does not display any SdH effect for either 
orientation. In contrast, sample B exhibits strong oscillations for H// 
as low as 2.5 T, with faint oscillations also emerging for > 6 T. 
However, data acquired with an in-plane field H//[010] do not show 
any oscillations up to 4 T. Symmetry dictates that the plane of a 2D FS 
in LaAlOa/SrTiOa must lie parallel to the interface; any such FS will 
therefore lack states with out-of-plane momenta and cannot exhibit 
any SdH effect for in-plane fields. It is therefore immediately clear 
that the oscillations which we observe with H//[110] must originate 
from an anisotropic 3D FS. 

For H//[110], the SdH oscillations in sample B are sufficiently 
pronounced for us to extract the effective mass and the Dingle 
temperature (a measure of the scattering) from their temperature- 
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Figure 1 | Shubnikov-de Haas oscillations and effective mass of the high 
mobility electron gas below the interface, (a), Oscillatory components of 
the MR RxxiH j_ J/)/ Rijg in samples A and B. The background i^^^was 
determined using a polynomial fit to the raw RxxiH) data for 1 .5 T < Hj^ < 
9 T and 2 T < H// < 4 T. Data for sample B were acquired at — 2.4 X 
10^^ cm"^, i.e. Vg = 350 V. The vertical scale for these graphs is identical, 
(b), Temperature variation of the in-plane SdH oscillation amplitude in 
sample B iVg= 350 V). (c), Amplitude suppression with increasing 
temperature for four in-plane SdH MR oscillation peaks/troughs Rp^i T) 
(indicated by arrows in (b)). Symbols represent oscillation amplitude data- 
points, while the four colour-coordinated lines are least-squares fits using 
equation 1 to calculate m*. Due to the low resistance of sample B (~1 Q), 
the SdH effect for H ± (001) creates voltage fluctuations of ~1 nV, which 
is the noise threshold in our measurements. Consequently, oscillations 
were only visible below T = 0.1 K and m* could not be determined for 
these carriers. 

dependent amplitude (Fig. lb). The magnitude of the oscillatory 
resistance is given by: 

(i^..(H,r)-i^,,(H,T))/4i^,,=exp( ^^)^:^/sinh(^^) (1) 



eHH 



eHH 



where R^g is the background resistance. Fitting this equation to the 
oscillation amplitude (Fig. Ic) yields = 1.24 ± 0.1 and To = 
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1.4 ± 0.4 K. is similar to values previously reported for the q 
LaAlOs/SrTiOs 2DEG^'^, although our is lower which implies a 
higher carrier mobility in our heterostructures. To estimate this ^ 
mobility, we initially calculate the Hall mobility fin = l/n2D^Rxxy 5 
where Rxx(yg = 0) = 0.28 Q/D and we assume single-band trans- ^ 
port. This yields an exceptionally high Hall mobility fin — ^ 
32000 cm^V~^s~\ setting a new record for pure LaAlOa/SrTiOa ^ 
and rivalling the best epitaxial SrTiOs films^\ o^^ 

In order to justify such a high mobility, we evaluate the Drude 
scattering time T^r = m^fiH/e = 23 ps, which is more than an order 
of magnitude greater than the Dingle scattering time = fil2%kBTu ^ 
= 0.87 ps. An alternative estimate of the scattering time in sample B 
may be extracted from the field at which a SdH effect first appears, 
using the quantum oscillation emergence condition cOc^sdh ~ 1 
(where cOc = Be/m* is the cyclotron frequency and B the magnetic 
field strength). For H//[110], oscillations are visible above 2.5 T: this 
corresponds to Tsdh = 2.8 ps, which is also shorter than t^^ sug- 
gested by our high ^h- It is likely that four factors contribute to this 
disparity: firstly, all scattering events suppress quantum oscillations 
and contribute to t^, while only back- scattering influences T^r and 
the Drude conductivity. Similar variance between To and T^r can be 
seen in other LaAlOa/SrTiOa heterostructures^. Secondly, the finite 
thickness of the conducting channel in our heterostructures may 
postpone the emergence of any SdH effect, until the applied field is 
sufficiently large for the diameter of the cyclotron orbits to fall below 
this thickness. Thirdly, superconducting fluctuations at fields below 
—2.5 T effectively "short-circuit" our heterostructures, reducing our 
ability to probe transport from carriers deeper below the interface. 
Finally, our single-band estimate for /Hh is an over- simplification, 
since multiband transport is expected for carrier densities above 
the Lifshitz transition in LaAlOs/SrTiOs^^'^^. A three-band approxi- 
mation to the field- dependent Hall coefficient (see Supplementary 
section 2) suggests a minority contribution from a high-mobility 
band with jih ~ 8000 cm^V"^s"\ The total number of conduction 
bands in our heterostructures and their field- dependent mobilities 
remain unknown, so we cannot obtain a more precise value for the 
mobility of these quantum- oscillating carriers. However, it is clear 
that our SdH effect, resistivity and Hall data all indicate the presence 
of a high-mobility band with an anisotropic FS and /Hh ~ 
10^ cm'V-'s-\ 

The fact that our measured is lower than than those reported 
for the LaAlOa/SrTiOa 2DEG^'^ also suggests that the band whose FS 
generates the in-plane oscillations lies within an extremely clean 
region of our heterostructures, far from the cation defects and mag- 
netic scattering expected at oxygen -deficient PLD-grown LaAlOs/ 
SrTiOs interfaces. To determine the location of these high-mobility 
carriers more precisely, we examine the evolution of the SdH oscil- 
lation frequencies with field-effect doping, obtained from the peaks 
in fast Fourier transforms (FFTs) of Rxx(H//,±) (Fig. 2a,b). The 
Onsager relation links the peak frequency F with the extremal area 

Sh 

S of the FS normal to the applied field via F = : since the size of 

the FS should be proportional to the carrier density, it is useful to 
compare F(Vg) with our experimentally- determined total as well 
as the superconducting critical temperature (which varies strongly 
with the local three-dimensional carrier density Hsd^^'^^). Once the 
interfacial carrier density exceeds n^D ~ 10^° cm"^, we expect a 
gradual suppression of SC leading to a dome in T^VgY'^; this is 
indeed observed (Fig. 2c). However, the in-plane oscillation fre- 
quency F// is independent of Vg, implying that the FS area S ± 
[110] responsible for these oscillations remains roughly constant 
upon field-effect doping. Furthermore, F//(Vg) displays no correla- 
tion with Tc(Vg) or n2D(yg)- the FS (and hence the density of states) of 
the SC band(s) is being influenced by field- effect doping, but the FS 
of the high-mobility band is not. Field-effect doping should have a 
similar effect on all occupied bands within the same spatial region. 
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Figure 2 | Evolution of the Shubnikov-de Haas oscillations and carrier 
density with increasing gate voltage, (a), Variation of the oscillating 
component of R^xiH//) for sample B (left panel) at Vg > 0, with FFTs of the 
raw data (right panel). For V^^ < 0 the noise level rises and it is not possible to 
identify oscillations: this is a well-known phenomenon and has been 
attributed to emergent inhomogeneit/. (b), Oscillating components of 
RxxiH±) for Vg > —100 V in sample B (left panel) with associated FFTs 
(right panel). The two peaks in the FFTs are indicated by grey and red 
arrows; for Vg = —100 V, the peaks merge, (c), Vg dependence of various 
properties of sample B, including Tc (above), SdH frequencies F//^2± and 
(below). Tc is measured from RxxiT) data (see Supplementary Fig. S2) and 
the errors in F// 2 ± correspond to the FFT peak widths at 80% of their 
maximum height (from (a),(b)). /72D(^g) obtained from the HaU coefficient 
foUows the values expected from the sample capacitance C{Vg) (see 
Supplementary Figs. Slb,S3b for raw capacitance and Hall data). We 
attribute the fall in n2D above V^^ = 350 V to charge-trapping deep within the 
SrTiOs. Tc(Vg) forms a dome: since d= 19 nm at Vg = 350 V and SrTiOs is 
SC for 5.5 X 10'" cm"' < < 5 X 10'° cm"' we estimate a maximum 
conducting channel thickness W ~ 20 fim due to the combination of carrier 
injection and electron gas decompression^. In practice, we anticipate 
W < 1 fim due to the extremely high at the interface which wiU locally 
suppress SC: even at Vg = 350 V, every carrier in sample B could be 
accommodated in merely 3 nm of SrTiOs doped at 0.5e"/unit cell. 

Therefore, the only possible explanation for this decoupling between 
TciVg) and F//{Vg) is that the SdH-oscillating electron gas must be 
spatially separated from superconductivity, i.e. the high-mobility 
carriers lie below the SC channel. 

The gate evolution of R^xiHj^ ) is very different from Rxx{H//), with 
two y^-dependent peaks appearing in the FFTs (Fig. 2b). One of 
these {Fi j_ , grey arrows) lies below 20 T and is suppressed for large 
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Figure 3 | Orbital occupancy and Fermiology calculations at high carrier densities in LaAlOs/SrTiOs. (a), Layer- dependent orbital populations for 
= 3 X 10^^, 3X10^^ and 8 X 10^^ cm~^. Data are plotted on a logarithmic scale with a lower cut-off of 10~^ electrons per Ti atom. (b),(c), Fermi surfaces of 
the interfacial d^y (red) and ^^^^2,72 bands (green) projected onto the (001) and (110) planes for = 8 X 10^^ cm"^ (thin lines) and 3 X 10^^ cm"^ 
(thick lines). The FS are calculated using a tight-binding model in which the parameters are fitted to bands calculated from first principles (see 
Supplementary Material for the band structure). For simplicity, our diagram ignores the hybridization which lifts the degeneracy at the band crossing 
points; taking this into account would split the doubly-degenerate FS into two. (d), Effective mass variation ml{kp) for ii^^z electrons within the (001) and 
(110) planes: thin and thick lines denote results at = 8 X 10^"^ and 3 X 10^^ cm"^. The planar angles are measured from the [100] and [001] axes 
respectively, ml (kp) for dyz electrons in the (001 ) plane is equivalent to ml (hp) for the d^z FS rotated by 90° around [001 ]. A sketch of the mass variation 
expected in the (110) plane for a bulk degenerate dxy,xz,yz is also shown (dashed brown line). (e),(f), Schematics illustrating the extremal FS orbits 
normal to magnetic fields along [001], [010] and [110] for a degenerate (i_x:y,xz,yz (e) and a (i^^,/^- dominated FS (f). Only H//[ 110] inthe dxz,yz-^ommsLted 
case probes the small, light FS cross-section whose presence we infer from our in-plane SdH oscillations. 



Vg. although this frequency seems too low to originate from the d^y 
interfacial bands (which form a larger FS at much lower n2D'^), spin- 
orbit splitting may create a series of small FS for high at the 
interface. The other peak (^2^, red arrows) mirrors n2D(V^) as Vg 
increases, saturating and broadening at —40 T for large Vg. This 
implies that F2 j_ also cannot arise from a d^^y 2DEG at the interface, 
since for back- gate doping at > 5 x 10^^ cm~^ the interfacial d^y 
occupancy should not change significantly: instead, carriers move 
deeper into the SrTiOs. It is therefore tempting to link this peak with 
the 40 T mode from de Haas-van Alphen experiments^^ on ^-doped 
bulk SrTiOs; however the light 3D band whose spherical FS was 
shown to be responsible for the 40 T oscillation^^ is only occupied 
for > 6.7 X 10^^ cm"^, by when SrTi03 already shows SC^^. Since 

< 20 nm for sample B^°, n^D ^ 5.5 X 10'^ cm"' below the SC 
channel, ruling out any occupancy of this light band. We therefore 
attribute F2 ± to the gradual population of the 3D FS from the first 
occupied band in bulk doped SrTiOs, which is formed by degenerate 
Ti 3dxy,xz,yz orbitals and remains approximately isotropic for such 
low n^D ^ 10^^ cm~^. 

It is clear that the in -plane SdH effect F// in our data is unrelated 
not only to F2±> but also to any previously reported 2D^'^ or 3D^^ 
quantum oscillations in LaAlOs/SrTiOs. Instead, our oscillations 
originate from a highly anisotropic FS (since there is no SdH effect 
for < 6 T), which occupies a clean intermediate region between 
the interface and the bulk. We estimate the minimum thickness of 

this region using the cyclotron radius Vg = : since a depth of at 



kF = y 27i;F/// (where is the magnetic flux quantum and we 

assume a spherical FS for simplicity), obtaining IVg ~ 140 nm at 2.5 
T. 

To understand the origin of these in-plane SdH oscillations, we 
calculate the evolution of the sub -interfacial orbital occupancy 
(which determines the FS symmetry) with increasing ^20- The 
majority of electronic structure calculations for LaAlOs/SrTiOs to 
date have only considered the first few layers below the interface for 
— 10^^ cm"^ and are of limited use in our heterostructures. We 
have therefore performed first-principles calculations of the depth- 
dependent band structure in LaAlOs/SrTiOs for n2D = 3 X 10^^, 3 X 
10^^ and 8 X 10^^ cm"^, specifically chosen to approach our experi- 
mental n2D in samples A, B (y^ ~ 0) and B ( > 0) respectively. Our 
calculated orbital occupancies are plotted in Fig. 3a and can also be 
seen in Fig. 4a-c: although computational power limits us to con- 
sidering the first 10 unit cells below the interface, this is already 
sufficient to reveal the FS anisotropy responsible for our in-plane 
SdH effect. 

The central result from these calculations is a crossover from d^y to 
dxz,yz occupancy as we move away from the interface. Close to the 
interface and for small ^20^ d^y states dominate due to quantum 
confinement, as expected^"'^^. The absence of a clear SdH signal from 
the 2D dxy interfacial FS in sample A is due to scattering from local 
moments^° and the large Rashba spin-orbit coupling; we note that 
there are no reports of a 2D SdH effect in FM LaAlOs/SrTiOs in the 
literature. The important new result from our calculations is the 
creation of a conducting "tail" deeper below the interface for large 



least 2r^ is necessary to establish SdH oscillations for H//(001), we use n2Dy with a disproportionate occupation of dxz,yz orbitals. For 
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example, the d^z,yz'-^xy ratio in layer 9 for /I2D = 8 X 10^^ cm"^ is 
2.8 : 1, significantly greater than the 2 : 1 expected in bulk SrTiOs. A 
recent study of top-gated SrTiOs also hints at a low density "tail" of 
carriers persisting over at least 50 Ti02 layers, independently of the 
total n2B^. While the majority of carriers occupy tightly-bound 
bands close to the interface, the back-gate field in our sample B 
should reduce the quantum confinement and expand the "tail" still 
further into the SrTiOa: this competition between confinement and 
decompression is responsible for the weak variation of 
(Fig. 2c). We therefore identify a dxz,yz- ^ommated FS as the source 
of our in-plane SdH effect. 

The strong asymmetry in our observed SdH effect (i.e. the absence 
of oscillations for small H^) may be explained by considering the FS 
geometry. In Fig. 3b,c, we plot the calculated (001) and (110) extremal 
cross-sections of the interfacial FS at /I2D = 3 X 10^^ and 8 X 
10^^ cm"^. The elliptical cross-section of the dxz,yz FS implies that 
our previously- calculated Vg will be scaled by kp[ooi]/kp[iio] = 0.73, 
reducing the minimum thickness over which the electronic structure 
deviates from that of bulk SrTiOa to —100 nm. Furthermore, the 
variation in \kp\ across the FS drives a corresponding modulation in 
the effective band mass m^, shown for the (001) and (110) planes in 
Fig. 3d. Electrons in the (001) plane are significantly heavier and hence 
more easily scattered: therefore, SdH oscillations will only emerge for 
Hi^H//. Our measured F// ~ 25 T is clearly too small to originate 
from the large interfacial FS projections in Fig. 3c: instead, our in- 
plane oscillations are generated by a similarly- shaped smaller FS dee- 
per below the interface, where is lower. The overall symmetry of 
the dxz,yz FS does not vary significantly with depth and hence our 
effective mass argument justifying the suppression of oscillations for 
Hj^ remains valid. In the (110) plane, the average band mass of 

the carriers is ml = 2(^ml^^^Q^ml^QQ^^^ j {m\^^^^^^ m\^^^^ =0.7 m„ 

which only allows for a small electron-phonon coupling k ~ 0.8 when 
compared with our measured = 1.24 (since m* = (1 + /l)mp. 
However, we note that SdH experiments on both LaAlOs/SrTiOs and 
n-type SrTiOs heterostructures have persistently yielded small effec- 
tive masses^'^'^^. 

Identifying the role of m\(kp) in determining the emergence of 
SdH oscillations allows us to make a profound statement regarding 
the shape of the in -plane oscillating FS. In Fig. 3d, we sketch the 
approximate m\ dependence in the (110) plane expected for a degen- 
erate (bulk-like) d^y^z^yz FS. Here, the m\ variation is similar to that in 
the (001) plane, though with a 180° rather than 90° period. We 
attribute the absence of oscillations for small to the presence 
of heavy carriers in the (001) plane: therefore, the emergence of 
oscillations at small H// implies that m\ cannot rise significantly at 
0°. Consequentially, the FS within this >100nm sub -interfacial 
region must be flattened along the [001] direction in comparison 
with the bulk, i.e. the d^y^^^^y^ degeneracy is lifted and the d^^^y^ orbi- 
tals are shifted to lower energy. To illustrate this point further, in 
Fig. 3e,f we sketch d^y^^z,yz and (i;^^^^^- dominated Fermi surfaces, com- 
paring the shapes of their extremal orbits perpendicular to [001], 
[010] and [110]. The low-frequency SdH oscillations which we 
observe with H//[110] must originate from a FS whose extremal 
orbits are composed exclusively of light carriers (i.e. the FS cross- 
sectional area must be small): it is clear that this condition is only 
satisfied for the (i;^^,^^- dominated FS. 

Discussion 

What is the physical origin of this change in the FS? We note that the 
shape of our d^^^y^ FS is similar to that calculated by Mattheiss^^ using 
a crystal-field parameter D which was subsequently shown to be too 
large^^. Since D is related to the tetragonal structure of SrTiOs, our 
renormalised electronic structure may result from strain effects at the 
interface - such as the compression from the LaAlOs layer - which are 
known to influence the 2DEG^^. Studies of the 2D-3D crossover in ^- 



doped SrTiOs films^^ (in which strain should be absent) have not 
revealed the dominated intermediate FS which we observe; 
nevertheless it remains unclear whether a long-range interface- 
induced change in D or the spin- orbit coupling is responsible for 
our results. Finally, our determination of the FS orbital character 
assumes the SrTiOs tetragonal c-axis lies parallel to [001]: since 
orthogonal tetragonal domains are expected for T < 105 K, this 
may not initially seem plausible. However, an offset surface potential 
exists between domains with c//[001] and c//[100] in LaAlOs/ 
SrTiOs, requiring substantial charge transfer to equalise the chemical 
potentiaP^. This increases the carrier density in domains with cH 
[001], so transport predominantly occurs within these regions. 
Previous transport studies of La-doped SrTiOs have also indicated 
a prevalence of [001] -oriented domains^^. 

We summarise the evolution of the LaAlOs/SrTiOs interface with 
in Fig. 4, where we schematically represent the spatial distri- 
bution of SC together with the approximate n^,^ variation and our 
calculated depth-dependent d^y and d^z^y^ orbital occupancies 
(Fig. 4a-c). At low carrier densities (Fig. 4a), d^y orbitals dominate 
and the charge is concentrated within a few unit cells of the interface. 
Electrons in the top Ti02 layer tend to locaUse^°'^\ creating an 
inhomogeneous patchwork of FM zones above a narrow 
(<20 nm) SC channeP". 

As ^20 increases (Fig. 4b), FM and SC both remain present at the 
interface. However, a high-mobility d^^^y^ "tail" of minimum thick- 
ness 100 nm develops below the interface, generating an anisotropic 
3D FS which exhibits SdH oscillations for small in-plane fields. 
Together, the appearance of this SdH effect, its independence from 
i^ioiVg) and TdVg), and its absence in small perpendicular fields 
indicate that dxz,yz orbital occupancy is favoured over d^y to a depth 
of at least 120 nm below the interface. Unfortunately, it is not pos- 
sible to accurately determine the maximum depth reached by this 
"tail", since the carrier density very close to the interface (where we 
expect the majority of the carriers to reside) is unknown. However, 
our data do enable us to comment on the vacancy penetration 
depth, which we already believe to be small since the capacitance of 
our B-type samples is comparable to values seen in annealed hetero- 
structures. The high electron mobility within the "tail" region is 
primarily a consequence of the low carrier density (which leads to 
a small FS and low effective mass), but a lack of crystal defects (e.g. 

vacancies) below the interface may also play an important role. 
Recently, ultra-high mobility carriers (fin ~ 50,000 cm^V~^s~^) have 
been observed in SrCu02-capped LaAlOs/SrTiOs heterostructures, 
in which vacancy formation is suppressed^^. This suggests that 
although the carriers in our B-type heterostructures originate from 
O^" vacancies, these vacancies may be confined close to the interface 
(or in the LaAlOs layer) while the electrons which they donate are 
redistributed deeper within the SrTiOs. This concept is supported by 
the absence of any parasitic SrTiOs surface conduction in our hetero- 
structures (whose presence would be expected in the case of deep O^" 
vacancy penetration), as well as theoretical work which indicates that 
O^" vacancies preferentially inhabit the LaAlOs surface rather than 
the interface^^ Ideally, future theoretical work should examine the 
evolution of the electronic structure in the "tail" as a function of O^" 
vacancy density and location. It also remains to be determined 
whether the absence of superconductivity from the "tail" region is 
merely due to a sub -critical carrier density, or if the dxz,yz orbital 
character also plays some role. 

At the maximum ^20 which we are able to simulate (Fig. 4c), only 
the top Ti02 monolayer at the interface still has a d^y character, with 
dxz,yz States dominating below. We illustrate the effects of a back-gate 
electric field in Fig. 4d: as Vg increases, the carrier density in the 
superconducting channel rises and a shift to the overdoped side of 
the superconducting dome occurs (as seen in Fig. 2c). In parallel, 
electrons in the "tail" decompress away from the interface due to 
band-bending from the electric field, migrating hundreds of nano- 
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a Sample A b Sample B, Vg ^0 c Sample B, Vg > 0 
n,n = 3x10i3cm-2 n,n = 3x10^4 cm-2 /?,„ = 8x101* citt^ 




Figure 4 | Electronic structure, orbital character and depth-dependent phase emergence at the LaAlOs/SrTiOs interface, (a-c), Schematics illustrating 
the phase and carrier distributions for sample A (a), B at low Vg (b) and B at large Vg (c). The approximate local carrier density within the SrTi03 is 
indicated by the colour shading: pale yellow denotes undoped bulk SrTiOs, while higher carrier density regions are either red {dxy)y green (d^z^yz) or 
brown (degenerate dxy,xz,yz) depending on the dominant orbital character, d^y ferromagnetism is present regardless of the carrier density, but remains 
tightly confined to the interface^°, i.e. at the top of the red zone. For comparison, we also plot the calculated orbital occupancies for the first ten Ti02 layers 
below the interface at = 3 X 10^^ (a), 3 X 10^^ (b) and 8X10^^ cm~^ (c). The index "1" refers to the Ti02 layer closest to the interface. Red and green 
bars correspond to d^y and dxz,yz orbitals respectively, (d). Qualitative illustration of the depth-dependent influence of field-effect doping on the local 
carrier density, which we deduce from our transport data. The line colour indicates the variation in the dominant orbital character with depth. 
(e),(f). Cartoon Fermi surfaces of the sub-interfacial high-mobility dxz,yz electron gas (e) and the bulk doped SrTiOa (f) which develops progressively for 
large n2D' 



metres into the bulk. This migration creates the 3D FS responsible for 
the SdH oscillations which we observe with H±(001), whose fre- 
quency scales with the total carrier density. Between the interface and 
the bulk, the carrier density of the (i^cz,^^- dominated region remains 
roughly constant: electrons which it "loses" to deeper-lying bulk 
states are replaced by electrons from the interface. The presence of 
a large carrier population below the interface results in a screening of 
the electric field, thus explaining the relatively small increase of d to 
19 nm at Vg = 350 V compared to (i > 40 nm reported at much 
smaller back-gate fields in the literature^^. Finally, Figs. 4e,f display 
exaggerated sketches illustrating the evolution of the FS as we move 
deeper into the SrTiOs, from dxz,yz domination (Fig. 4e) to a gradual 
recovery of dxy,xz,yz degeneracy (Fig. 4f) over a lengthscale > 120 nm. 
While the microscopic origins of this long-distance evolution are still 
unclear, our work shows that functional oxide devices can reliably 
hope to profit from a renormalised electronic structure tens of nano- 
metres away from a symmetry-breaking interface. 

Methods 

Two series of LaAlOs/SrTiOs heterostmctures, "A" and "B", were grown using a 
standard pulsed laser deposition system manufactured by Twente Solid State 
Technology B.V., equipped with a reflection high-energy electron diffraction 
(RHEED) facility. We use 0.5 mm thick commercial 5X5 mm SrTiOg (001) "STEP" 
substrates from Shinkosha: these are HF-treated for Ti02 termination and cleaned by 
the manufacturer, then vacuum-packed for shipping. We do not perform any addi- 
tional surface cleaning or annealing prior to deposition: the substrates are loaded 



directly into our PLD chamber, which is subsequently evacuated to base pressure 
( < 10 ~^ mbar) prior to back-filling with 10"^ mbar O2. The substrate is then heated to 
growth temperature (800°C). Series A and B both feature 10 unit cells of LaAlOa, 
deposited using a total incident laser energy of 9 mj focussed onto a 6 mm^ rect- 
angular spot. The O2 pressure and substrate temperature were maintained at 10"^ 
mbar and 800°C respectively for both sample series throughout the deposition pro- 
cess. Subsequently, A-type samples underwent an annealing stage: after cooling to 
500°C at 10"^ mbar, the O2 pressure was increased to 0.1 bar. The temperature was 
held at 500°C for 30 minutes before natural cooling to 20°C in 0.1 bar O2. In contrast, 
B-type samples were cooled naturally to 20°C in 10"^ mbar O2. 

To fabricate Hall bars on these LaAlOs/SrTiOs films, we first defined contact pad 
areas using photolithography with AZ52 14 photoresist. 2 nm Ti followed by 8 nm Au 
were evaporated directly onto the LaAlOa surface; the remaining photoresist was 
removed by soaking in acetone for 30 minutes, then rinsed in IPA. Sample B also had 
an Au-Ti back gate deposited across the entire base of the SrTiOs substrate prior to 
fabrication. The Hall bars were defined using a similar photolithography process and 
the Hall bar mesas etched using a dry Ar ion technique (at a slow rate of 1 A s"^ to 
avoid any substrate heating). The Hall bar width was 80 [im and the voltage contact 
separation 660 \im. Multiple Hall bars were fabricated on each 5X5 mm substrate: 
tests showed that the Hall bars were electrically isolated from each other (thus ruling 
out any parasitic conduction from the SrTiOs surface) and displayed similar transport 
properties (indicating that our heterostructures are homogeneous). Prior to mea- 
surement, the Hall bars were mounted in thermally- conductive chip-carriers, with 
electrical contacts made using 10 [im Au wires ball-bonded to the Au-Ti contact pads. 

Transport data were acquired in a cryogen-free dilution refrigerator, using an AC 
technique with two digital lock- in amplifiers and a current source outputting 500 nA 
at 19 Hz. This value was chosen to maximise the signal-to-noise ratio whilst mini- 
mising sample heating below 0.1 K. Our noise threshold is approximately 1 nV. The 
substrate capacitance was measured with femtoFarad sensitivity for gate voltages up 
to 500 V using a General Radio 1621 manual capacitance bridge. All results presented 
in this work were qualitatively reproducible over a 6 -month period comprising 
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numerous cool-downs of both samples. A total of 6 "A-type" and 4 "B-type" het- 
erostructures were fabricated in our laboratory using identical "recipes" to those 
detailed above: all samples displayed similar behaviour to those discussed in the 
present work. 
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